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We discuss the relation between continuum bound states (CBS) localized on a defect, and surface 
states of a finite periodic system. We use the transfer matrix method to model an experiment of 
Capasso, and find all continuum bound and anti-bound states. We compute the rate for intra- 
subband transitions from the ground state to the CBS and derive a sum rule. Finally we show how 
to improve the confinement of a CBS while keeping the energy fixed. 
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I. INTRODUCTION 

In 1992, i^n experiment of F. Capasso et al. 
demonstrateda the existence of well localized continuum 
bound states (CBS) in a semiconductor superlattice con- 
sisting of one thick quantum well surrounded on both 
sides by several GalnAs-AlInAs well/barrier layers con- 
structed to act as A|/4 Bragg reflectors. As suggested by 
Lenz and Salzmarn, the central well was made double 
the width of the lattice wells, to 3jCti clS 3t A/2 Fabry- 
Perot resonator. Subsequently, Webem studied the exis- 
tence of such states using the transfer matrix method. 
Among other things, he showed that the Bragg con- 
dition need not be very well, satisfied for a confined 
state to exist. B. Sung et al.a have also studied above 
threshold confined states, in a different material system 
GaAs/AlGaAs. Radovanovic et al.El have written several 
papers describing mainly numerical methods for tailor- 
ing a heterostructure to produce CBS with equal energy 
spacing, to facilitate second—harmonic generation. Fi- 
nally, Xue-Hua Wang, et al.B have discussed the parity 
sequence of sub-threshold bound states localized on a de- 
fect, and the transition rates between them, but neglect- 
ing the variable effective mass. 

In this paper we provide further insight into the phe- 
nomenon of CBS by relating them to surface states, 
whose existence was elucidated by Shockley in a de- 
servedly famous papeiiH. An infinite periodic system, il- 
lustrated by line (c) of Fig. ^ allows Bloch states with 
the periodicity of the lattice. If the system is truncated 
on one side, or on both sides, (line (b)) then one can 
discuss scattering states with energies above threshold, 
and bound states below threshold. The transfer matrix 
method is well adapted to discuss such a periodic system. 
For convenience we will use the notation of our previous 
papersQcl. 

Among the bound states of a finite periodic array are 
the Bloch states whose amplitude is spread more or less 
uniformly over the lattice, and the surface states whose 
density is confined to the ends. The former usually oc- 



cur in the allowed energy bands of the infinite lattice, 
while the surface states necessarily occur in the forbidden 
bands. Their wavefunctions decay exponentially outside 
the array and like ±e~^ from cell to cell inside, where 9 
is the imaginary part of the Bloch phase. 

Another way to truncate an infinite lattice is to cut 
it in the middle and pull the two halves apart. This 
introduces a defect (line (d) of Fig. |l|) which may be 
either a well or a barrier. As emphasized by WeberEI, 
the condition for a state localized in the central gap of 
an infinite lattice is that the boundary condition at the 
edge of the defect matches to a decaying Bloch eigenstate 
of the unit cell: that is, the wave function will decay by 
±e^^ from cell to cell. There are also anti-bound states, 
where the match is to the growing eigenstate. In either 
case, the Bloch phase must be complex, so such states 
exist only in the forbidden Bloch zones (FZ). 

In section 2, we set up the (generalised) transfer ma- 
trix for a system with position-dependent effective mass. 
In section 3 we apply it to the Capasso experiment. We 
determine the width of the central well/barrier to pro- 
vide a CBS at a desired energy. In sec. 4 we discuss 
the relation between these states and states in a box, 
illustrated in frames (a) and (f) of Fig. |l]. In Sec. 5 
we show how to compute the transition rate from the 
ground state to continuum states in the neighbourhood 
of the CBS, and derive a sum rule which implies that the 
total strength depends very little on the number of cells 
involved. Finally, in Sec. 6 we propose new experimental 
configurations that improve the confinement of a CBS 
while keeping the energy fixed. 



II. TRANSFER MATRIX 

Consider solutions of the Schrodinger equation with a 
variable effective mass 
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At a discontinuity of the potential or of the effective mass, 
both ip'(x)/m{x) are continuous. Let v{x) and 

w{x) be any two independent solutions at the same E. 
Then the modified Wronskian 
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takes a constant value. We choose C = 1 by suitable 
normalization. 

For convenience we will discuss the situation where 
the entire system has reflection symmetry. Then it is 
sufficient to consider solutions only for x > and reflect 
them in the origin. In many cases, including the specific 
examples discussed by Weber, Capasso and Wang et al., 
it is possible also to choose the unit cell to have reflection 
symmetry about its midpoint. We will do this when it 
aids in the analysis. 

We factorise m(x) — mem*{x) into the bare mass rUe 
and the dimensionless m*. The constant /2m = 3.81 
eV A" converts from energy to length units. In Weber's 
model both the potential and the effective mass are piece- 
wise constant functions, but the method is valid even if 
V{x) and m*{x) vary continuously within the cell. 

Consider a solution v{x) which has value 1, slope at 
the left edge of the unit cell, and another solution w{x) 
which has value and slope m*(0) there. The transfer 
matrix for the unit cell of the lattice is then 
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Here the solutions v,w without argument are evaluated 
at the cell edge, x = d — 0, and tp is an arbitrary solu- 
tion. In a periodic system, W{d) depends only on the 
length of the cell, not its position. Since ip and -p' jm* 
are continuous at interfaces, to move one interval further 
to the right one simply multiplies again by the appropri- 
ate transfer matrix. Any discontinuity in the derivative 
is automatically taken into account. 

Eigenvalues of the transfer matrix satisfy 



A — 2 cos -f 1 = where 
2cos0 = 1xW — V + w'/m* 



(4) 



and when the Bloch phase is a real angle, they are A = 
g±z0 Generally, raising the energy from the potential 
minimum, one is in a forbidden band where | cos0| > 1. 
In this region of energy, (j) = i9 is imaginary. Following 
this, the first allowed band occurs within which (j) in- 
creases from zero to tt. Then in the next forbidden band, 
with cos(/) < —1, (j) — n + i9 becomes complex. In the 
p'th forbidden band, (j) = pi: + iO and the eigenvalues are 
A = (-)Pe±^ 



A. Surface States 

Because we have assumed that parity is a good quan- 
tum number, states of the whole system will have either 
even or odd parity. Suppose that the infinite array is 
truncated so that there are N cells to right of the ori- 
gin as in Fig. 1, frame (b). Then the condition for a 
bound state is that the wave function at the right edge of 
the array matches to a decaying solution outside. (Here 
we suppose constant potential outside, but that can be 
changed trivially.) 
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For an even bound state, ^ — v{x), on the right hand 
side acts on (1,0) while for an odd state w{x) it 
acts on (0, 1). This gives the log-derivative 
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where s = 1 (2) for even (odd) states, and m*{d) is the 
effective mass inside the edge of the last cell, while m%^^^ 
is the value outside, and E = Vout — / {2mm*^^). By 
construction, the M^-matrix is real, so the energy E must 
lie below the external potential Vout ■ On the other hand, 
a surface state can exist only when is in a forbidden 
zone, with complex <j) = pn + i9. In such a zonecl, 

W{Nd) = W'^id) 



sinh^ 



sinh ( 
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Eqs. ^, ^ allow one to search for energies where surface 
states occur. 



B. Continuum Bound States 

Suppose that the infinite periodic array is cut at the 
origin and an extra well of width 2c is placed between 
the two sections, as in frame (d) of Fig. |. Let T(c) 
be the transfer matrix (as in eq. |[) that takes the wave 
function from the origin to c. Its columns are the even 
and odd parity solutions within the central well. In order 
for a CBS to exist, the first (or second) column of T(c) 
must match to a decaying eigenstate of the unit cell of 
the semi-infinite array to the right. In other words, one of 
the columps of T{c) must satisfy the eigenvalue equations 
for W{dp: 
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Either of these equations can be used to search for CBS. 
If they are satisfied with real A being the smaUer (larger) 
eigenvalue, then a CBS (or an anti-bound state, (ABS)) 
exists at that energy. The only difference between them 
and surface states is the numerical value of the boundary 
condition that has to be satisfied where the lattice meets 
the defect, or the surface. 



III. CAPASSO- WEBER EXAMPLE 

The array constructed by Capasso can be modelled as 
a sequence of potential wells of width a — 16A, depth 
= and barriers of width b = 39A, height VJ, = 500 
meV. The energy-dependent effective mass in each layer 
is given byO 



= 0.043(1 + {E- V^)/E^] 
= 0.073(1 + {E- Vb)/Eb) 



(10) 



where E^ = 0.88 eV, and Eb = 1.49 eV are the effec- 
tive band gaps of the InGaAs wells and AlInAs barrier 
materials. 

In this example, the potential is piece-wise constant, 
so the transfer matrix can be constructed from factors of 
the type 



T(c) 



cos kc 
— fc sin kc/m* 



m* sin kc/k 
cos kc 



(11) 



where fc^ — 2mgm*{E — Vc)/h'^ is the wave number inside 
the layer and to* is the effective mass there. In contrast 
to Weber, we find it convenient to define the unit cell as 
a half- well of width a/2 on each side of a barrier of width 
b: then 



W{d) = T{a/2)T{b)T{a/2) 



(12) 



is the transfer matrix of a symmetric unit cell, and has 
both diagonal elements equal to cos (j). Panel (a) of Fig. 
^ shows a part of the array. 

In Fig. ^, we show the eigenvalues in the forbidden 
zones. To guide the eye, they are connected by straight 
lines across the allowed zones. There is an even index 
forbidden zone up to 307 meV in which the eigenvalues 
vary enormously. The odd FZ from 387 to 641 meV shows 
a smaller but still appreciable range of variation. This is 
reduced still further in the third FZ from 881 to 960 meV. 
The smallest values of \e~^\ occur near the middle of the 
forbidden zones. 

The transfer matrix from the origin to c + Nd is 
{d)T{c). For any non-zero c the central well con- 
stitutes a defect in the superlattice. If ^ oo, the 



argument of Sec. 2B applies, and the wave function 
will only be localized near the origin when one of the 
columns of T{c) is an eigenvector of W{d) with the de- 
caying eigenvalue ±e~^. The left hand side of eq. ^ 
(cf. eq. [ll] ) is either —k^tankwc/m^ (for an even 
state), or kw cotkyjc/mlj (odd state) . But the identity 
cot(a;-|-7r/2) = — tan(a;) means that the solutions for odd 
states can be found simply by adding 7r/2 to the value 
of fciuC of an even state solution. Given any solution, an- 
other one which differs only by the number of nodes in 
the central well, can be obtained by adding ir/k^ to c. 
Hence, it is the differences in width of the central well 
that go by half- wavelengths, not the whole width. 
By solving 



Wii{d) — A TO^ 

Wl2{d) 'ho 

W2iid) 

W22{d)-\ k^ 
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we obtain the values of c in Fig. Q where we plot the 
solutions from the principal branch of the arctan func- 
tion. In the first (even) FZ, lines 2 and 4 reach the c = 
axis at the beginning of the allowed zone, at 641 meV. In 
our convention, c — makes the central well have total 
width a — 16A, the same as in the infinite SL, so one 
has a periodic system in which no confined state exists. 
That is why the loci of CBS or ABS terminate on the 
line c = at either end of a forbidden zone. 

The horizontal dotted line at c = 8A picks out the odd 
CBS state at 563 meV found by Weber and Capasso. In 
this case, the central well is double the width of the wells 
in the lattice. Very close by at 577 meV is an even ABS 
(line 3). In the next odd-parity FZ there is an odd parity 
CBS at about 890 meV. 

Fig. H allows one to read off the energies of the CBS 
(ABS) for any choice of width c. Line 2 is just line 1 
augmented by 7r/(2fcu,); similarly for lines 4 and 3. Neg- 
ative values of c down to —a/2 are meaningful, as they 
simply shrink the width of the central well down to zero. 
The straight line segments across the allowed zones are 
drawn just to connect related solutions. 

Because the choice of a flat central potential leads to 
the tangent function in eq. |l^, one can invert to obtain c 
directly. For a more complex potential one would have to 
search for solutions of eq. |^, but nothing would change in 
principle. Any potential that produces the same bound- 
ary condition at c gives the same solution, with only the 
details of the wave function between ±c changing. For 
example, Wang et al.H used a two-segment well/barrier 
combination as their central well. Its boundary value still 
varies between ±oo, so one must expect similar solutions, 
albeit at slightly different energies. _ 

As explained by ShockleyQ, (see alsoEl) a unit cell which 
is more attractive at the edges than in the centre, is just 
the type which favours the occurrence of surface states. 
Compared to a well of constant average depth, the ener- 
gies of even states are raised and those of odd states, low- 
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ered. This leads to band mixing in the forbidden zones, 
making it possible for the decaying eigenstate of the cell 
to match to a decaying exponential with a positive log 
derivative of tp at the outer edges. In the Capasso device, 
it is the odd state in the central well which has the posi- 
tive log derivative where it matches to the lattice on the 
right, making the CBS. Viewed from the lattice, at the 
left edge the match is to a fast decreasing wave function 
in the central well, which must vanish at the origin. 

To illustrate the close relationship between CBS and 
surface states, in Fig. || we show in the upper panel 
the right side of a six cell lattice with the central well 
of width 2a, having the CBS at 563 meV. In the lower 
part we show a three cell lattice with a surface state at 
the left edge decaying into the potential barrier, chosen 
so the slope is exactly the same for both wave functions. 
The only difference is that one state goes to zero at the 
origin to make an odd parity wave function, while the 
other decays exponentially; within the lattice they are 
identical. 

As a second example, we replaced the central well by 
a central barrier, (see Fig. j^b), similar to the system 
proposed by Lenz and Salzmancl. In this case the unit 
cell is taken to be a barrier of width 6/2 on either side 
of a well of width a. The diagonal elements of W{d) 
are interchanged but its trace is unchanged, so the band 
structure is identical. In Fig. ^ we show the correspond- 
ing results for the extra barrier-width c. As before, c ~ 
corresponds to a periodic system and no CBS is possi- 
ble, as this point comes at the edge of an allowed band. 
The central barrier pushes up the energies of even parity 
states, so now an even-parity CBS is found in the neigh- 
bourhood of 550 meV, with c ~ 20A. There is a nearby 
odd-parity ABS if c > llX. 

To see whether the thicknesses chosen by Capasso were 
optimal, we varied b and a. For fixed well width a, the 
area enclosed by the eigenvalues of the odd FZ in Fig. ^ 
expands continually as b is increased, even to 200A. The 
lowest absolute value of the smaller eigenvalue |A_| gets 
smaller and smaller: see Table |. However, this simply 
means that the wells become progressively decoupled as 
the distance between them gets larger. Similarly, for fixed 
barrier width b = 39 A, |A_| decreases steadily as the well 
width increases, while the energy of the CBS moves down 
towards threshold. At a = 20A it becomes a true bound 
state below threshold. For a = 19A, |A_| = 0.2508 which 
is almost twice the rate of decay of probability (which 
depends on |Ap) from cell to cell within the lattice, as 
compared to a = 16A. The energy of the CBS is then 
at 506 meV. The values k^a = 1.810 and hb = 0.426 
are even further from the Bragg condition than for the 
original parameter set. This confirms Weber's finding 
that a high degree of confinement does not require exactly 
satisfying the Bragg condition. 

For a central barrier, the behaviour with increasing 
barrier width is similar to what was found above for a 
central well: see Table ||. Here it is an even parity CBS 
that falls in energy from 778 to 503 meV as b increases 



from 20 to 100 A. But again, this is a result of decoupling 
of the wells. For fixed b = 39A, there is a shallow mini- 
mum of |A_| = 0.351 at a = 17A and the energy of the 
CBS is at 574 meV. Again, k^a = 1.757 and kbb = 1.505 
are far off the Bragg condition. 

However, the conclusion is very different if both a and 
b are varied while keeping the energy of the CBS fixed. 
An example is shown in Table III. The minimum value 
of I A- 1 is obtained with both ky^a and ktb — 7r/2. It 
should be expected because, once we fix the energy of 
the state, the effective mass is also fixed. Then the op- 
timization of k^a and k^b proceeds exactly as for an 
energy-independent Kronig-Penney potential, for which 
the Bragg condition is optimal, as one can easily show 
analytically. 



IV. RELATION TO STATES IN A BOX 

Kalotas and LeeEl considered the states obtained by 
enclosing a finite number of cells between infinite walls. 
(See frames (a) and (f) of Fig. |l|for an illustration.) This 
discretizes the continuum, so all states become discrete. 
Well localized states that decay quickly enough will be 
scarcely affected by the walls. States spread over the 
whole lattice will become a discrete set maintaining a 
similar character. An ABS whose magnitude grows away 
from the origin will be squeezed against the walls of the 
box. 

In order for the finite N system to agree with Weber's, 
we add an extra half-well or barrier at the right edge of 
the array. Fig. |^ shows an example where we have taken 
iV = 10 cells on each side of the origin. This is to be 
compared with the CBS at E = 563 meV of the c = 8A 
example of Capasso and Weber. Even with just three 
cells on either side, the state is hardly shifted from its 
position in the infinite array. 

In Fig. ^ we show the spectrum of box states as a 
function of N, again for the c = 8A central well case. 
The energies of the single-cell states change little as more 
cells are added. The new states that appear fill up the 
allowed bands. To understand this, it is convenient to 
consider a system with a hard wall at the origin, then 
N identical cells, followed by a hard wall at the right. 
The allowed wave functions are those that vanish at the 
origin (odd-parity states of the symmetric system), and 
the hard-wall boundary condition at the right edge is 
■0(x = Nd) = 0. In view of eq. ||, this requires that the 
element Wi2{Nd) = 0. Since in an allowed band 
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This shows that bound states can occur in either of two 
ways. First, as single-cell bound states, where the second 
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factor vanishes. These states have A'' nodes, and the wave 
function vanishes at every cell boundary. Alternatively, 
the combinatorial factor sinA^0/sin<^ might vanish, and 
in an allowed band there are N — 1 such states where 
N(j) = TOTT, m = 1,2, ...N — 1. The single-cell state may 
occur in a forbidden zone, but the others can only occur 
for real 0, in an allowed band. For non-zero c one has to 
multiply (d) from the right by the additional transfer 
matrix T(c) for the central well, and then the simple 
factorization won't be exact. In practice the states tend 
to remain in the allowed band all the same. 

Fig. 1^ provides another example of the similarity of 
CBS and surface states. In the upper panel (a), the CBS 
lies in the middle of the first of the first forbidden zone, 
while in the lower panel there are two such states, one 
derived from allowed band a and the other from band /3. 



V. TRANSITION RATES 

Introducing the vector potential into the Hamiltonian 
eq. |l| leads to the excitation operator 



eA 1 1 



—iheA 
2m„c 



■S 



(16) 



with dimensions of energy. In defining the operator 5, 
(dimensions inverse length) we have factored out the bare 
electron mass, leaving only the dimensionless effective 
mass (m* ~ 0.06) inside. The vector potential A is as- 
sumed to be a function of x, so it commutes with the 
mass. By invoking the Coulomb gauge we make it com- 
mute with the momentum as well. 

According to the Golden rule, the transition rate is 



27r / eAh 



h \2meC 



<^f\S\^.>fpiEf), (17) 



where p{Ef) is the density of final states. The factors be- 
fore the matrix element have dimensions length squared 
energy per second, and these are omitted from our cal- 
culations. The matrix element squared times the density 
of states is therefore (energy length-squared)"^, and this 
is what we plot in Fig. 10 for example. After integrating 
over energy, we have units for the total strength. 

In the Capasso experiment, the ground state has even 
parity, so its derivative is odd, and transitions are allowed 
only to odd parity excited states. Also the lattice is finite 
rather than infinite, so the transitions are to states in 
the continuum. In the neighbourhood of the CBS, the 
continuum wavefunction has a large normalization inside 
the central well, and this causes the transition rate to 
peak at or near this energy. 

We again consider the case of a central well of half- 
width c surrounded by N cells on each side. But for 
convenience, in this section we follow Weber by defining 
W{d) — T{a)T{b) as having a barrier of width b on the 
left followed by a well of width a. Then N will be the 



number of wells additional to the central well. The odd 
parity excited state with Dirac delta-function normaliza- 
tion has wave function 

'i/ f{x) ^ Bosinik^^x) , |a;| < c 

= ^sm{kbX + 5) x>c + Nd. (18) 



where £;-yfc = n^kl/{2 mml) measures the energy above 
the top of the barrier in the asymptotic zone, and fc^ = 
2mm*^E/h^ is the wave number inside the central well. 
Using the transfer matrix W{d) to cross N cells gives 

1 f sm{kb{c + Nd) + 5) 
y/^ \ cos(fc{,(c -f- Nd) -f 5) 

sin k,,,c 



l^in COS knjC 



(19) 



where u,^ = k^/m*^, Vb = kb/nib and ly(^) = W^{d), 
so that: 
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ta.TL{kb{c + Nd) + 5) = 

w[f'' sin fc^c -f ly^'^Vtu cos k^c 
W2i^ sin k^c + W22^i'w cos fc^c 



(20) 



determines d and 

^sin(fcb(c-|-iVd) + S) 

AN) 



Bq(Wi" ' sinfc^c -I- W^2^iyw cosku,c) (21) 

gives the normalization Bo- Note that the matrix el- 
ements of W^^'^ can be easily computed from those of 
W{d) using eq. ^ (or eq. |l^ when the Bloch phase 4> is 
real). One need not solve explicitly for the phase shift 
5{E) because only |i3oP is required to compute the tran- 
sition rate, and the identity sin^ z = tan^ z/{l + tan^ z) 
can be used in eq. 

With the above equations we can construct the wave- 
function f{x) as follows. Wave functions v{x), ^/{x) in 
a unit cell of the lattice are defined in cq. ^. Within the 
r'th cell following x = c, (r = 1, 2, ...) the wave function 
is 

"^rix) = Arv(x — c — rd + d) + Brw{x — c — rd + d) 

(22) 

and from the matching at x = c we have: Ai = 
Bosin/c^jC, Bi = Bt)Vyj / VbCOskyjC. In general. 



Ar+l 

Br+i/m* 



W{d) 



Br/m* 



l,...(iV-l). 

(23) 



This gives values of "^rix) in each cell and allows 

the calculation of the matrix element in eq. |l^. 
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The ground state wave function is computed in a sim- 
ilar manner. Inside the central well it is 



^'o(a;) = No cos(fc^x) 



X < c 



(24) 



and at a; = c + Nd it must match to a decaying expo- 
nential, as in eq. If the state is well bound it is a 
good approximation to use the Kalotas state which van- 
ishes at the edge of the lattice or the Weber state that in 
principle extends to infinity. The normalization constant 
A^o niust be computed by summing the normalization in- 
tegrals from every cell as well as from the central well. 
If one integrates < v\v >, < v\w >, and < w\w > over 
the unit cell, then it is just a matter of multiplying these 
integrals by the coefficients in the r'th cell and summing. 

Because the effective mass depends both on position 
and momentum, it is not obvious how to evaluate the 
matrix elements of the transition operator S. It is rea- 
sonable, in the term p/m* to let the p act on the excited 
state ijjE and interpret the effective mass as being at that 
energy. Conversely, in the {l/m*)p term, where the p acts 
on the intial state, we use the ground state effective mass 



TOq- Then, 



< iIje\S\iPo > 



dipt 



1 



dx m*{E) 
' m*{EQ) dx 



Tpodx ■ 



dx 



(25) 



Since tp' /m* is continuous at interfaces between wells and 
barriers, the integrand is continuous, despite the jumps in 
m* . When the effective mass m* is piece- wise constant, 
we can evaluate the integral over a series of intervals of 
constant m* (here interpreted as the m*{x,E)). 

The squared matrix element, including the density of 
states factor, is plotted as a function of energy in Fig. 
[lOL which is to be compared with Fig. 2 of Sirtori et 
aOJ. (However, their figure has normalized the peak in- 
tensity to unity in each case, obscuring the fact that the 
integrated strength is constant.) As the number of side 
wells increases, the computed excitation function rapidly 
becomes very narrow. (The experimental width is much 
larger than theory, indicating that something else is hap- 
pening.) It shows that even a small number of cells is 
sufficient to give a well confined state. We also find in- 
creasing strength in the second allowed band near 700 
meV as cells are added. The integrated strength under 
the main peak (from 500 to 640 meV) varies only a few 
percent. 

To illustrate how the continuum wave functions evolve 
in the region of the CBS, we show in Fig. ^ four cases 
spanning the energy range. It can be seen that as one 
passes over the CBS energy at 563 meV, an additional 
node appears in the wave function. Away from this reso- 
nance, (panels (a), (b)) the wave function consists mainly 
of the growing solution in the lattice, so the amplitude 
is largest at the outside edge (where it is fixed, accord- 
ing to eq. fsh. Close to the CBS, (panel (c)), there is 



a large component of the decaying solution, making the 
amplitude in the central well large. Increasing the energy 
again, (d) , brings back more of the growing solution. At 
577 meV, the position of the ABS, only the growing solu- 
tion would contribute. If we had more than iV = 3 cells, 
the effects would be even more pronounced. 

While the peak in the transition strength becomes very 
narrow as the number of lattice cells N increases from 
zero to 3, the integrated strength is almost constant. 
This can be understood from the sum rule which follows 
from eq. and is discussed in the appendix. The total 



A2. In addition to the 



strength, M2, is defined in eq. 
integral over the continuum, when there are N Bragg re- 
flectors on either side of the central defect, there will be 
N discrete odd-parity bound states, which must also be 
included in the sum. Typically these account for some- 
thing like 6% of the total strength. These odd bound 
states are shown in Fig. |l2|, for the case N = 2. In this 
figure, the wave functions are remarkably similar inside 
the region where the ground state is large, so they give 
almost equal contributions to the sum rule. 

Turning now to the results, one has to distinguish be- 
tween the no-reflector case and the TV-reflector case. In 
the former, the integrated transition strength (ITS), M2, 
is about lAA^^. The strength is very broadly distributed 
above threshold; see line (a) in Fig. ^ We have looked 
at = 1 to 5 cells on each side. For these cases, the ITS 
is around 1.35A^^, of which about 0.09 comes from the 
bound odd-parity excited states. As N increases, the ITS 
fluctuates only a little. The strength remains highly con- 
centrated into the CBS peak (about 80%), the remainder 
being spread quite widely (6% in the bound states and 
14% in the continuum). 

Because the CBS peak becomes so narrow, we estimate 
the integral under it by assuming a Breit-Wigner shape, 
and deducing the height and width from the calculations. 

The sum rule, calculated according to eq. AS, is 
always about 6 to 10 % higher than the ITS, if we 
(arbitrarily) set the doorway state mass to be at 
E = EcBS = 563 meV. The main term in eq. A3 is 
proportional to (l/mg -I- l/m^)^, so we can easily ad- 
just the doorway energy to ensure that the sum rule will 
agree with the ITS. We call this the effective doorway 
energy Ed- With no reflectors Ed is about 740 meV 
(535 above the ground state energy). This is reasonable 
since the excitation strength is very broadly distributed 
above threshold. With one reflector (on each side), Ed 
drops to 660 meV; but then it slowly rises, at least up to 
N = 5, where it reaches 670 meV. 



VI. OPTIMAL CBS CONFINEMENT 

In this section we discuss some general principles for 
designing a well confined CBS. Consider a general unit 
cell of width d, within which the potential and the effec- 
tive mass are arbitrary functions of x. (Since we will fix 
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the CBS energy, this allows for energy dependence of the 
effective mass.) Now let us arbitrarily divide the cell into 
two parts so that the widths a and b add to d, and denote 
by W^, the transfer matrices of the two parts. When 
a is on the left, we have ^W^W, with elements 



(26) 



If the whole array is symmetric about the origin, there 
will be two type-a portions together at the origin, and the 
sequence of potentials is ...baba...ba\ab...abab.... When we 
look for odd parity confined states of such an array, it 
is equivalent to putting a hard wall at the origin, and 
solving only the right side. The wave function at the 
edge of the first cell will be, using eq. || 



■0 /m* 



<2 



X 



(27) 



The second equality holds if we imagine an infinite array, 
and demand an eigenstate with the wave function in each 
cell differing only by a factor A. This wave function will 
vanish at both x = and d, and in the second cell, the 
value of tp' /m* will differ by a factor A from the first. The 
condition that must be satisfied is = 0, and then 
A = W22- Writing these out in terms of the component 
transfer matrices gives 



<2 = W^,W^2 



W^2W^2 = 



Mi 



W^^W^2 



and 



ly,", = A 



'22 



'2 1 



2 2 



Mi 



"^12 



= {WIM2 
WI2 



- W^22<i) 
Wa{a) 

Wb{b) 



W?2 



Wt2 



(28) 



In the last step we have used the form of W as in eq. ^. 

What this tells us is that to make the eigenvalue as 
small as possible, we must make the Wa{x) solution as 
small as possible at the right edge of the a part-cell, and 
conversely, Wb{x) as large as possible at the right edge of 
the b part-cell. 

The above is true for any division of the cell into two 
parts. In the system studied by Capasso et al., the logical 
division is into the two layers of GalnAs and AlInAs. In 
that situation, the off-diagonal elements have the form 
= — s\n.kcc/V(, where = kc/m* (c = a, b) is the 
velocity. If the Bragg reflection condition holds, then 
sinfcc = ±1 and A = —v^jva is just the ratio of velocities 
in the two parts of the cell. (This is analogous to the 
problem of waves on a string, with one part thin and the 



other thick. At the join, the displacement y(x) is contin- 
uous, and the ratio of the slopes y'{L)/y'[R) is the ratio 
of the velocities squared.) One sees that in the a-cell, 
the solution Wa{x) rises to the value Wa(a) — l/t'a, while 
starting from d and moving backwards through the bar- 
rier region, the corresponding solution falls to the value 
Wb{—b) — —Wb{b) = —X/vb- Normalizing the b solution 
to ensure continuity at x = a requires the factor A. 

The result eq. ^ is quite amazing because in the gen- 
eral situation where the potential and effective mass vary 
arbitrarily, the dividing line can be placed wherever you 
wish. To make the eigenvalue small, one must make the 
"odd-parity" solution w{x) as large as possible every- 
where in the second part-cell and as small as possible in 
the first. As observed by Weber, the first aspect can be 
achieved by choosing an energy just above the barrier 
(small kb). To meet the Bragg condition, this forces a 
large 6, and the linear variation of Wb{x) over the bar- 
rier leads to a large wave function at a; = 6. That is 
why lowering the energy of the CBS in general improves 
confinement. 

However, our aim is to improve confinement while keep- 
ing the energy of the CBS fixed. To do so, we are guided 
by Shockley's explanAtion for the existence of surface 
states, as discussed iiJa. We split the a part-cell into two 
sections, ai, a2, making the left side more attractive, and 
the right side less so {Vwi < Vu, < ¥^,2)- Then the greater 
curvature of Wa near the origin, balanced by less curva- 
ture to the right, will lead to a smaller value of Wa{a), 
even if the average attraction is the same. We leave the 
6-cell fixed, but a similar strategy with less repulsion on 
the right side, can obviously be employed there. 

To illustrate this, we have selected a set of heterostruc- 
tures based on quaternary alloys GuxIrii-xAsyPi-y, 
lattice-matched to InP {x = 0A68y). We took infor- 
mation from Figs. 1.17 (for band gaps) and 1.20 (for 
effective masses) of Swamiiia|than and MacranderliJ. For 
the band offsets, S. Adachillj gives 



AEc = 268y + 3?/2 
AEy = 502y - 152?/2 



(29) 



The baad alignments are also discussed on p. 87 of 
DaviesE3. Putting together this information, we arrived 
at the following set of parameters: 



y=1.0 m; = 0.043 = 0.880 eV 



y 
y 



0.5 
0.0 



0.061 
0.081 



Eg 

En 



1.080 eV 
1.360 eV 



(30) 



where_Eg are effective band gaps in the sense of Nelson 
et aie. 

In this way we can have conduction band po- 
tential steps of 125 meV (from Ga. 47 /n 53 As to 
Ga. 23/^.77 As. 5P5 or 250 meV to InP (y = 0). These 
are a quarter (denoted Q) or a half (denoted H) step up 
to the 500 meV barrier of ylZ.48/n.52As. 
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As our baseline, (denoted {Q,Q) below), we take the 
a-well to consist entirely of Q (y= 0.5) material, so the 
potential floor at 125 meV is 375 meV below the barrier. 
The barrier width was fixed at 44. 3A which satisfies the 
Bragg condition. For a width a — 15.82A, the CBS is at 
63 meV above the top of the barrier, as in the original 
experiment. The eigenvalue A = —0.4026, which is not as 
favorable as in the original work because the well is not so 
deep. (Capasso et al. selected the materials to have the 
greatest possible well-barrier potential difference.) The 
potential properties are summarized in Table IV, top line. 



Next we divide the a-well into two parts, one of 
GalnAs {y = 1) and the other of InP. We adjusted 
the widths a i an d 02 to keep Eqbs fixed. In the second 
line of Tabic denoted (0, H)i, the two part wells are 
of equal width 9.353A, and the eigenvalue is A = —0.396. 
In the third line, denoted {0,H)2, the deeper well has 
width ai ~ 10.5QA, and the shallower part a2 = 7.91lA, 
giving A — —0.374. This may seem a small gain, but we 
shall see that the improvement is significant. 

We then computed the CBS properties for a finite ar- 
ray based on the above materials, with results shown in 
Tables 0, ^ and Fig. |l3|. The transition strength is 
significantly narrower and more strongly peaked for the 
split- well examples. The strength to the bound states in 
the split-well cases was only 25-30% of that of the (Q,Q) 
reference case, and the strength to the continuum states 
was much larger. The results are summarized in Table 

where the number of cells means the number of Bragg 
reflectors placed on each side of the central defect. 

The differences in the total strength are reflected in 
the portion concentrated in the CBS peak. The decay 
constant strongly influences both the sharpness of, and 
the area under, the peak. The cases with a lower (in 
magnitude) decay constant have narrower, and larger (in 
terms of area) peaks in the transition strength curve as 
seen in Table ^ 

Overall, the decay constant has a signiflcant effect on 
both the total strength, and its continuum and bound 
state contributions individually. A lower value of the 
decay constant (in magnitude) results in better confine- 
ment of the CBS as is evidenced by the width of the peak 
in the transition strength curve. We conclude that the 
split-well strategy can produce much better confinement 
of the CBS. It should be feasible to confirm this method 
of improving the confinement of CBS, experimentally. 



VII. CONCLUSION 

We have shown that Continuum Bound States are 
closely related to surface states, because both arise as a 
result of perturbing an infinite periodic system. The re- 
sults of Webeitl concerning the experiment of Capasso et 
al. were verified. However, while we confirm that a CBS 
can exist even when the Bragg conditions are not well 
satisfied, we note that for this type of potential, one can 



prove analytically that optimal confinement is achieved 
by the Bragg conditions. In addition, by enclosing a fi- 
nite array in a box, we have traced the evolution of the 
Bloch states in the allowed bands as the number of Bragg 
layers is increased. 

We have derived a sum rule, within the conduction 
band only model, that explains the integrated transition 
strength from the ground state to the continuum in the 
region where the CBS exists. About 70% of the sum 
rule strength is concentrated in this region. Finally, we 
have identified the factors that allow one to improve the 
confinement of a continuum bound state, and proposed 
a way of testing this. 
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APPENDIX A: SUM RULE FOR 
INTRA-SUBBAND TRANSITIONS 

Sirtori et al.lll discussed the sum rule for excitations to 
the CBS within a two-band Kane model. Here we limit 
our discussion to what can be done within a conduction- 
band only model, eq. |l|. The difficulty which arises in this 
context is that having an energy-dependent effective mass 
means that the Hamiltonian does not have a complete 
orthonormal set of excited states, so the sum rule can 
only be approximate. 

We take the operator to be 

- ihS = {p— + —p) (Al) 

TTlE mo 

(cf. eq. |l6[) Then the sum rule is 

h'^M2 ^n^ < o\s''s\Q > 

= < 0\{p — + p){p + — p)|0 > 

mo ms niE niQ 

= < 0\p— — —p\Q > + < Q\p— — —p\Q > - 

mo mo mo ms 

in < 0|p— ( — ))'|0 > + < 0\p— — —p\0 > + 
Too mE mE mo 

ih < 0|( — )'— p|0 > + < 0\p— —p\0 > + 

niE mo mE mE 

ih < 0|( — )'— p|0 > -ih < 0|p— ( — )'|0 > 
mE mE mE mE 

+?i2<0|(— )'( — )>> , (A2) 

mE mE 
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Ma = 



' + —rdx+ < 0|( — )'( — )'|0 > 

mo niE niE ruE 



Mo 



^ 1 

mo 



rriE niE 



dx 



(A3) 



To obtain this expression we moved the p-operators until 
they act on the ground state wave function directly. In 
the case of a constant effective mass, too = m^; = 1, only 
the first integral survives. In this case the sum rule must 
be exact, and we found close agreement between the sum 
rule expression, and direct integration: 

/dk 
< ^o\S^\^Pe > ^ < ^£1-51^0 > dE 

where 



dk 1 2m , 2{E-Vb). 
dE = 2kl^'^"'^^ + ^^> 



(A4) 



which includes non-parabolicity. The wave number k is 
defined by the energy above the barrier. At high energy 
the density of states tends to a constant, rather than 
going to zero, as it would for constant mass. The case 
without energy-dependence can be recovered if i?b — > oo. 
The density of states factor was tested by doing the in- 
tegrals either over energy or over wave number k. 

When we introduce x-dependence to the effective mass, 
the terms involving the derivative of (l/m^) contribute. 
Because in this model the mass is piecewise constant, the 
derivative is a Dirac delta function times the discontinu- 
ity in (1/to^). The integral in the last line of eq. A3 is 
then a sum of values evaluated at the layer edges. 

When we introduce energy-dependence as well, both 
here, and in the first integral, factors such as V'o/'^b 
discontinuous, because the mass TOe is taken at one en- 
ergy and the ground state wave function at another. To 
resolve this ambiguity we took the average of the two 
values on either side of the discontinuity. For these ma- 
terials, the well and barrier masses are similar, so it is 
not a large uncertainty. This is the stage at which the 
sum rule can only be approximate. Moreover, we need a 
prescription for the energy E at which we evaluate uie- 
Thinking in terms of the doorway state approximation, 
initially we took the CBS e nerg y. 

The second term in eq. A3 involves the square of a 
Dirac delta function, and is undefined. We simply omit 
this contribution. 



J. Radovanovic, D. Indjin, V. Milanovic and Z. Ikonic, Solid 
State Comm. 110 (1999) 339-343. 

Xue-Hua Wang, Ben- Yuan Gu, Guo-Zhen Yang and Jian 

Wang, Phys. Rev. B58 (1998) 4629-35. 

Wm. Shockley, Phys. Rev. 59 (1939) 319-323. 

D.W.L. Sprung, Hua Wu and J. Martorell, Am. J. Phys. 

61 (1993) 1118-24. 

D. W. L. Sprung, J.D. Sigetich, Hua Wu and J. Martorell, 
Am. J. Phys. 68 (2000) 715-22. 

D.F. Nelson, R.C. Miller and D.A. Kleinman, Phys. Rev. 
B 35 (1987) 7770-3. 

T.M. Kalotas and A.R.Lee, Euro. J. Phys. 16 (1995) 119- 
23. 

C. Sirtori, F. Capasso, J. Faist, and S. Scandolo, Phys. 
Rev. B 50 (1994) 8663-74. 

V. Swaminathan and A.T. Macrander, "Materials aspects 
of GaAs and InP based structures", Prentice Hall, New 
Jersey, (1991) pp. 1-42 

S. Adachi, "Physical Properties of III-V Semiconductor 
Compounds", Wiley Interscience, New York (1992) pp. 96- 
9. 

John H. Davies, "Physics of Low-dimensional Semiconduc- 
tors", Cambridge U.P. (1998). 



TABLE L CBS properties varying widths a and b sepa- 
rately, for a central well. 



Central Well 



a 


b(A) 


E (meV) 


A- 


kiij a 


kbb 


16 


16 


686 


-0.565 


1.876 


1.011 


16 


20 


645 


-0.499 


1.797 


1.105 


16 


30 


589 


-0.389 


1.684 


1.275 


16 


39 


563 


-0.327 


1.632 


1.383 


16 


50 


545 


-0.277 


1.595 


1.483 


16 


60 


534 


-0.244 


1.574 


1.557 


16 


75 


525 


-0.210 


1.555 


1.647 


16 


100 


516 


-0.173 


1.537 


1.767 


16 


125 


512 


-0.150 


1.528 


1.862 


16 


150 


509 


-0.135 


1.522 


1.942 


16 


200 


505 


-0.115 


1.516 


2.072 


10 


39 


666 


-0.601 


1.149 


2.320 


11 


39 


651 


-0.545 


1.243 


2.202 


12 


39 


635 


-0.493 


1.332 


2.071 


13 


39 


618 


-0.444 


1.416 


1.926 


14 


39 


600 


-0.400 


1.494 


1.765 


15 


39 


581 


-0.361 


1.566 


1.585 


17 


39 


544 


-0.298 


1.693 


1.148 


18 


39 


525 


-0.273 


1.750 


0.860 


19 


39 


506 


-0.251 


1.801 


0.426 



^ F. Capasso, C. Sirtori, J. Faist, D.L. Sivco, S- N.G. Chu 

and A.Y. Cho, Nature 358 (1992) 565-7; C. Sirtori et al. 

App. Phys. Lett. 61 (1992) 898-900. 
^ G. Lenz and J Salzman, App. Phys. Lett. 56 (1990) 871-4. 
3 T.A. Weber, Solid State Comm. 90 (1994) 713-6. 

B. Sung, H.C. Chui, E.L. Martinet and J.S. Harris Jr., 

Appl. Phys. Lett. 68 (1996) 2720-22. 
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TABLE II. CBS properties varying widths a and b sepa- 
rately for a central barrier. 



Central Barrier 



a 


b (A) 


E (meV) 


A- 




kbb 


16 


16 


778 


-0.626 


2.058 


1.273 


16 


20 


715 


-0.555 


1.934 


1.372 


16 


30 


620 


-0.430 


1.747 


1.496 


16 


39 


577 


-0.354 


1.661 


1.539 


16 


50 


549 


-0.289 


1.605 


1.561 


16 


60 


535 


-0.246 


1.575 


1.570 


16 


75 


523 


-0.200 


1.550 


1.575 


16 


100 


513 


-0.153 


1.530 


1.577 


16 


125 


508 


-0.123 


1.521 


1.577 


16 


150 


506 


-0.103 


1.516 


1.576 


16 


200 


503 


-0.077 


1.511 


1.575 


10 


39 


602 


-0.443 


1.070 


1.785 


11 


39 


597 


-0.416 


1.169 


1.734 


12 


39 


592 


-0.396 


1.268 


1.689 


13 


39 


588 


-0.379 


1.367 


1.648 


14 


39 


584 


-0.368 


1.465 


1.609 


15 


39 


581 


-0.359 


1.563 


1.573 


17 


39 


574 


-0.351 


1.757 


1.505 


18 


39 


571 


-0.352 


1.854 


1.472 


19 


39 


568 


-0.356 


1.950 


1.439 



TABLE III. CBS properties for fixed energy 560 eV. 



Energy fixed at 560 meV 



a 


b(A) 


Energy 


A- 




kbb 


10 


77.29 


560.00 


-0.591 


1.016 


2.671 


13 


65.14 


560.00 


-0.393 


1.321 


2.251 


14 


58.25 


559.99 


-0.345 


1.423 


2.013 


15 


49.70 


559.98 


-0.318 


1.524 


1.717 


15.4 


45.98 


560.00 


-0.315 


1.565 


1.589 


15.45 


45.51 


560.01 


-0.315 


1.570 


1.573 


15.455 


45.465 


560.01 


-0.315 


1.5706 


1.5713 


15.457 


45.449 


560.00 


-0.315 


1.57084 


1.57065 


15.46 


45.42 


560.00 


-0.315 


1.5712 


1.5597 


15.5 


45.05 


560.00 


-0.315 


1.575 


1.557 


16 


40.45 


559.96 


-0.319 


1.626 


1.397 


17 


31.99 


560.01 


-0.348 


1.727 


1.106 


18 


25.25 


560.01 


-0.398 


1.829 


0.873 


19 


20.16 


560.02 


-0.460 


1.931 


0.697 


21) 


l().:5:i 


5(i().();i 


-().-",27 


2.():!:i 


l)..")64 



TABLE IV. Optimized three-layer potentials indicating 
the changes in the widths ai, a2, w{a), Va — fca/m*, 
and eigenvalue A. In all cases b = 44.2851450A, 
Wb{b) = 2.1450935. and Ecbs = 563.0 meV. 



case 


(11 


(1-2 


"■(") 




A 


{Q,Q) 


7.9109 


7.9109 


0.8636 


1.1579 


-0.4026 


(0,i/)i 


9.3535 


9.3535 


0.8495 


0.9092 


-0.3960 


{0,H)2 


10.4995 


7.9109 


0.8031 


0.9092 


-0.3744 
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TABLE V. Evolution of the transition strength {A ^) with increasing number of Bragg reflectors. 





Bound 


Continuum 


Total Strength | 


# Cells 


(Q,Q) 


iO,H)i 


(0,//)2 


(Q,Q) 




iO,H)2 


(Q,Q) 


{0,H)i 


(0,//)2 














0.912 


1.279 


1.324 


0.912 


1.279 




1.324 


1 


0.108 


0.0267 


0.0319 


0.762 


1.188 


1.234 


0.870 


1.215 




1.266 


2 


0.117 


0.0281 


0.0334 


0.759 


1.188 


1.233 


0.876 


1.216 




1.267 


3 


0.118 


0.0282 


0.0335 


0.754 


1.183 


1.229 


0.872 


1.211 




1.263 


4 


0.118 


0.0282 


0.0335 


0.746 


1.176 


1.224 


0.864 


1.204 




1.257 


5 


0.118 


0.0282 


0.0335 


0.745 


1.174 


1.225 


0.863 


1.202 




1.258 



TABLE VI. Total strength under the CBS peak: dependence on number of Bragg reflectors. 





Peak Height eV^A-^ 


Width r (eV) 


Peak/Total Strength (%) | 


Cells 


(Q,Q) 


(o,/i-)i 


{0,H)2 


(Q,Q) 


(0,H)i 


{0,H)2 


(Q,Q) 


(0,H)i 


{0,H)2 


1 


27.5 


38.2 


45.6 


1.73 * 10"^ 


1.49* 10"^ 




1.33 * 10"^ 


85.8 


73.7 


75.5 


2 


165. 


247. 


328. 


2.42 * 10"* 


2.07 * 10"'^ 




1.68 * 10"* 


71.6 


66.1 


68.4 


3 


1026. 


1609. 


2383. 


3.81 * 10"* 


3.16 * 10"" 




2.33 * 10"" 


70.5 


66.0 


69.1 


4 


6300 


10240 


16970. 


6.18 * 10"^ 


4.92 * 10"® 




3.24 * 10"* 


70.8 


65.7 


68.7 


5 


38950 


65450 


121300 


1.00 * 10"^ 


7.76 * 10"® 




4.54 * 10"® 


71.0 


66.4 


68.8 
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b) 

c) 




e) 



1 1 

I _ _ _ _ _ 



FIG. 1. Schematic drawing of an infinite array (line c), truncated to a finite array (b) and enclosed in walls (a); or with a 
defect (d), also truncated (e) and enclosed (f). 
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FIG. 2. A unit cell of the infinite lattice for A) a well as the central defect, and B) a barrier as the central defect. 
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FIG. 3. Eigenvalues vs. energy for a central well, c = 8A, showing the smaller magnitude eigenvalue A- (line 1) and the 
larger A+ (line 2). The straight lines across the allowed zones (where the A± are complex) connect related solutions. 
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FIG. 4. Central well width c versus energy showing the even CBS (line 1), the odd CBS (line 2), the even ABS (line 3), and 
the odd ABS (line 4). See also figure 6. The straight lines across the allowed zones connect related solutions. 
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FIG. 5. Central well width c versus energy showing a) the CBS and b) the ABS. Odd numbered lines are for even states, 
and vice versa. The straight lines across the allowed zone connect related solutions. 
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FIG. 6. Wave functions of a CBS (a) and a surface state (b) compared. 
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FIG. 7. Same as Fig. 4 but for a central barrier. (The straight lines across the allowed zone connect related solutions.) 
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FIG. 8. States in a box, for c=8A, also showing a) the CBS at 563 meV, b) the ABS at 577 meV, c) the CBS at 891 meV 
and d) the ABS at 946 mcV. At right, a labels the first allowed band from 307^387 rncV, d the second allowed band from 
641— >881 meV, 7 the third allowed band from 960— >1357 meV, and 6 the fourth allowed band from 1357 meV. Dotted lines 
show the potential cells, in eV; wave functions are dimensionless and are drawn with base line at the energy eigenvalue. 
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FIG. 9. Even (+) and odd (O) state energies (in a box) versus number A'" of cells to the right of the central defect. Panel a) 
is for a central well, and b) a central barrier. Bands labelled as in Fig. 8. 
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FIG. 10. The transition strength times density of states (units ev ^) for a) a central well of width 32A, and a central 
well surrounded by b) one cell, c) two cells, and d) three cells on each side. 
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FIG. 11. Four representative continuum wave functions, for c=8A, at energies E = 521 meV (a); 545 meV (b); 565 meV (c); 
and 573 mcV (d). Dotted lines show the potential cells, in eV; wave functions are dimensionless and are drawn with base line 
at the corresponding energy. 
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FIG. 12. The derivative of the ground state wave function (over m*) a) and the first and second odd-parity excited wave 
functions, b) and c), for a 5 well potential (two identical cells on each side of a central well), enclosed in a box, illustrating the 
similar overlap near the origin. 
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FIG. 13. Transition strength to continuum (units ev ^) for (a): case (Q,Q), (b): (0, and (c): (0, i?)2, showing the 
evolution of the CBS peak. Note change of scale between the upper panel for = 1 reflector and the lower panel for two. 
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